Measurement of the Transverse Polarization of Electrons 
Emitted in Free Neutron Decay 
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The final analysis of the experiment determining both components of the transverse polarization 
of electrons (oTi , crT2 ) emitted in the /3-decay of polarized, free neutrons is presented. The T-odd, 
P-odd correlation coefficient quantifying cjt2 , perpendicular to the neutron polarization and electron 
momentum, was found to be i? = 0.004±0.012±0.005. This value is consistent with time reversal 
invariance, and significantly improves both earlier result and limits on the relative strength of imag- 
inary scalar couplings in the weak interaction. The value obtained for the correlation coefficient 
associated with ar^, N — 0.067±0.011±0.004, agrees with the Standard Model expectation, pro- 
viding an important sensitivity test of the experimental setup. The present result sets constraints 
on the imaginary part of scalar and tensor couplings in weak interaction. Implications for parame- 
ters of the leptoquark exchange model and minimal supersymmetric model (MSSM) with _R-parity 
violation are discussed. 

PACS numbers: 24.80.-|-y, 23.40.Bw, 24.70.+S, ll.30.Er 
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I. INTRODUCTION 

In the Standard Model (SM) the description of free 
neutron decay involves only three parameters: (i) the rel- 
ative strength of axial and vector couplings, A — gA/gv, 

(ii) the first element of the quark mixing matrix, Vud, and 

(iii) a time reversal violating phase $. The much larger 
number of observables which became accessible in novel 
experiments at new generation neutron sources allows not 
only to contribute to the determination of those param- 
eters, but also to address some basic problems reaching 
beyond the SM. 

One of these is the incomplete knowledge of the physics 
of CP violation (combined charge conjugation and par- 
ity symmetry). The SM with the Cabbibo-Kobayashi- 
Maskawa (CKM) mixing scheme jlj accounts for CP vio- 
lation discovered in kaon and B-meson 0, 0] systems. 
It fails by many orders of magnitude to account for the 
most striking evidence of CP violation: the dominance 
of baryonic matter in the present universe Q . 

Supported by firm theoretical considerations and 
strong experimental evidence 0, the combined CPT 
symmetry is regarded to be a strict symmetry of na- 
ture. With this assumption CP violation is equivalent 
to time reversal symmetry violation (TRV) and as a re- 
sult is linked to microscopic reversibility and the princi- 
ple of detailed balance. No compelling evidence of TRV 
has been observed in experiments testing this principle 
in different nuclear reactions Q , and the T-violating am- 



plitude was found to be at most 10 ^ of the dominant 
strong interaction amplitude. 

To the most precise tests of time reversal invariance 
belong searches for electric dipole moments of elementary 
particles, atoms and molecules. Despite their impressive 
accuracy one obtains only upper bounds (2.9 x 10~^^, 
5.9 X 10~^^, 3.1 X 10~^^ e cm for neutrons [§], electrons 
[lOj and ^^^Hg atoms [ll|, respectively) which are still 
orders of magnitude away from the S M p redictions (e.g. 

_ e cm for the neutron Oll^), leaving a 

lot of room for new physics searches. 

The situation is more complicated in high energy ex- 
periments or in systems with heavy quarks contents. Also 
here it is possible to construct observables sensitive to 
TRV, but the sizable contributions of heavy quarks makes 
it difficult to disentangle between new physics and the 
SM induced effects. This was the case in the first direct 
observation of TRV in the kaon system in the CPLEAR 
[T3 | and KTeV experiments 19]. Only very recently the 
DO collaboration reported the observation of a charge 
asymmetry like-sign dimuon production in proton an- 
tiproton collisions at 1.96 TeV center of mass energy, 
which contradicts the SM at 3.2 standard deviations 16]. 

More than 50 years ago it was recognized that TRV 
may be tested also in various correlations accessible in 
nuclear or particles decays [13, [13 ■ Many systems have 
been investigated in this way including mesons [l^, lep- 
tons [lOj, baryons, and nuclei (see [2l| for a review). The 
measurement of the R coefficient in *Li decay, quanti- 
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fying the correlation between the spin of the decaying 
nucleus, the electron momentum and the electron spin, 
provides the most stringent direct limit on the imaginary 
part of tensor coupling constants of the weak interac- 
tion [22]. The discovery of new CP- or T- violating phe- 
nomena, especially in systems built of quarks of the first 
generation, with vanishingly small contributions from the 
CKM matrix induced mechanism, would be an important 
milestone. 

Free neutron decay plays a particular role in nuclear 
beta decay experiments searching for TRV. Due to its 
simplicity it is free from model dependent corrections as- 
sociated with nuclear and atomic structure. Further, fi- 
nal state interaction effects, which can mimic T violation, 
are small in this case and can, in addition, be calculated 
with a relative precision better than 1% [231. From a va- 
riety of correlation coefficients which may be built from 
vectors accessible in neutron decay, up to now only two 
have been addressed experimentally. First was the an- 
gular correlation between the neutron spin, the electron 
momentum and the neutrino momentum, referred to as 
D coefhcient in the literature. It is sensitive to the rel- 
ative strength {gA/gv) and phase angle {'^va) between 
axial and vector currents in weak interaction, and has 
been measured in several experiments 24-27]. At the 
current precision it provides the best limits to certain 
time reversal violating parameters appearing in standard 
model extensions with leptoquarks exchange, associated 
with a non-zero value of sin (PvA ■ 

In this paper we present the final analysis of the first 
measurement of another time reversal violating correla- 
tion coefficient in neutron decay - the R coefficient, and 
of the time reversal conserving N correlation, both asso- 
ciated with a correlation between the neutron spin, the 
electron momentum and its polarization. Being sensitive 
to the real and imaginary parts of scalar and tensor cou- 
plings of the weak interaction they provide information 
complementary to the D coefficient. 

This is the final report of the nTRV experiment com- 
prising data collected between 2004 and 2007. It super- 
sedes our previous result, presenting the methods used 
in the data analysis in more details and introducing a 
new i?-evaluation approach based on a "double" ratio 
method. The significant improvement in the accuracy of 
the determination of this correlation coefficient as com- 
pared to the result presented in [28*] is a consequence of 
two major extensions in the analysis of the existing data: 
(i) the analysis of an additional event class with backscat- 
tered electrons trajectories contained within the vertical 
plane, and (ii) improved determination of the effective 
analyzing powers of the applied Mott scatterers. Minor 
changes in the value of the P-even, T-even N correlation 
coefficient are the result of new effective analyzing pow- 
ers and the analysis of another event class which has also 
not been included in the previous analysis. 

The outline of this paper is as follows. In Sec.|TT]we in- 
troduce the measured observables and present the strat- 
egy of the experiment. The dedicated neutron beam line. 



detector setup and performance as well as the hardware 
trigger are discussed in Sec. Hm Sec. HV] presents the 
data analysis, discussion of systematic effects and de- 
scribes applied consistency checks. The results obtained 
and their implications on some extensions of the SM are 
compared with existing experimental data in Sec.|Vl and 
finally, conclusions are given in Sec. I VII 

II. CORRELATIONS IN NEUTRON ^-DECAY 

The electron distribution function for an experiment, 
in which the decaying neutrons are oriented, and electron 
energy, momentum ^E, p) and polarization are measured, 
is proportional to [l7|: 

M^((J),^,i?,p)a i + fe| + £. (^A^ + G^) + 

where m is the electron mass, J and (J) are the neutron 
spin and its expectation value, respectively, o- is a unit 
vector onto which the electron spin is projected and A is 
the /3-decay asymmetry parameter. Higher order terms in 
p, J and (T are neglected. N and R are the correlation co- 
efficients associated with uti and gt^ , respectively, where 
(7^2 is the transverse component of the electron polariza- 
tion perpendicular to the decay plane spanned by the 
neutron spin and the electron momentum, and axi is the 
component contained within this plane (Fig. [1} . As the 




FIG. 1; (Color online) Schematic representation of neutron 
decay. The decay plane containing the neutron polarization 
J, the electron momentum pe and the transverse component 
of the electron polarization ctti is indicated. 

terms involving G and Q are proportional to the longitu- 
dinal component of electron polarization which was not 
accessible in the presented experiment Eq. ([1]) reduces to: 

W{3,a,E,p)<x 1 + 6- 

+ P-(4+iV^ + i?^), (2) 

where P = (J) /J represents the average neutron beam 
polarization. It has been pointed out in [2l| that the 
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Fierz interference term (bm/E) affects most of the cor- 
relation measurements using neutron spin asymmetry to 
extract correlation coefhcients in a way, that the mea- 
sured quantity becomes 



X = X/{1 + (bm/E)) 



(3) 



with X — a, A, B, D . . . and the averaging is performed 
over the observed /3-spectrum. 

A. Final state interaction and exotic couplings 

Following [2^ the general Lorenz invariant interaction 
Hamiltonian density of nuclear beta decay can be written 
as: 

+ I /2 {pax^n){e(T\^{CT + C'j^lh)v) 
- {Plt^lbn){e-it,-i5{CA^C'Alb)v) 

(p75n)(e75(Cp-f C^75)i^) + H.c, (4) 

where Ci and C[ represent 10, in general complex cou- 
pling constants, which determine the symmetry proper- 
ties of the weak interaction. In the minimal formulation 
of the SM only vector and axial- vector interactions are 
present (Cy = Cy — 1, Ca — C'j^ — X) and all other cou- 
plings vanish {d = C,' = 0; i = S,T,P). With these 
assumptions, both R and N vanish at the lowest order in 
neutron decay, but acquire finite values when final state 
interactions are included: 



^FSI 
^^FSI 



A, 

E ' 
m 

-Qffs— - A, 
P 



(5) 
(6) 



where ais is the fine structure constant. For the energy 
distribution observed in the present experiment one ob- 
tains A^FSi ~ 0.068 and i?Fsi ~ 0.0006. This means that 
the i?FSi is tsLT below the sensitivity of this experiment, 
while the finite value of TV should easily by measured. 

Allowing for a small admixture of exotic couplings and 
keeping only terms linear in these couplings, one finds 
[17,] that: 



N - Npsi 

R — i?FSI 



-0.218 •3?(S') 
-0.218 •3(5') 



0.335 •3fi(r), (7) 
0.335 •!3(r), (8) 



where 



S = {Cs + C's)/Cv, 
T ^{Ct + C^)/Ca, 



(9) 
(10) 



are the relative strengths of scalar and tensor interactions 
with respect to the dominant vector and axial- vector cou- 
plings, respectively. Within these assumptions the coef- 
ficient b can be expressed as fvj^: 



and would affect the measured correlations following 
Eq. However, the additional terms are of second 

order in the contributions of exotic couplings and can 
thus be neglected. 

A non-zero value of the R correlation in neutron decay 
would signal the existence of a non-vanishing contribu- 
tion from imaginary couplings in the weak interaction, a 
new source of the TRV and, as a consequence, physics 
beyond the SM. 

The N correlation depends on the real part of the same 
linear combination of scalar and tensor couplings as R. 
However, the discovery potential of its measurement is 
strongly suppressed by a significant contribution of un- 
certainties connected with the evaluation of the final state 
interaction. 

There exist very few measurements of N and R corre- 
lations in general [l^ |33| , and only two in nuclear beta 
decays [H, |3l|. 



III. 



EXPERIMENT 



0.170 •5R(S') -I- 0.830 •3fi(r), 



(11) 



The key feature of the nTRV experiment is the abil- 
ity to measure energies and to track over relatively long 
distances electrons from neutron decay. This allowed for 
efficient use of one of the world strongest polarized, cold 
neutron beams as a source of electrons from neutron de- 
cay and for application of efficient electron polarimetry 
based on Mott-backscattering [32| . An additional advan- 
tage of this principle is the unique signature of relatively 
rare Mott scattering events which made them easily dis- 
tinguishable off- but also on-line from an overwhelming 
background of electrons accompanying this, very strong, 
neutron beam. 

Though a similar concept of electron detection has al- 
ready been applied in the measurement of the neutron 
life time and the A correlation coefficient [H, H^l, the 
present experimental setup outperforms the former ones, 
providing much more accurate reconstruction of both the 
electron trajectories and their energies. 

The experiment was performed at the SINQ facility of 
the Paul Scherrer Institute, Villigen, Switzerland. 



A. Cold neutron beam 

A dedicated cold neutron beam line has been con- 
structed for the present experiment at channel 51 of 
SINQ, leading directly to the cold moderator container. 
The container was filled with about 20 liters of liquid deu- 
terium at 25 K. Cold neutrons from the moderator were 
polarized in a 1.6 m long multichannel bender-polarizer 
[35| . and subsequently transported to the experimental 
area via a rectangular channel, farther referred to as con- 
denser (Fig. (2). Its convergent, vertical walls matched 
the 80 X 150 mm^ entrance beam cross section with the 
about 5.5 m distant 40 x 150 mm^ exit. The condenser's 
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FIG. 2: Schematic top view (not to scale) of the polarized 
cold neutron beam line arrangement. Cold neutrons from the 
liquid deuterium moderator enter from the left. 



main role was to increase the neutron density at the ex- 
periment, and to separate this area from a large back- 
ground of fast neutrons and gammas produced in the 
SINQ interior and in the polarizer. To minimize back- 
ground and neutron losses due to interactions with gas, 
the polarizer and condenser were enclosed in a vacuum 
chamber with 180 /im and 125 /im thick zirconium en- 
trance and exit windows, respectively. 

The application of carefully chosen, different kinds of 
supermirrors [sg'-'ss'I, with a critical reflection angle up 
to 3.3 times larger than that of natural nickel, applied 
in the polarizer and covering the walls of the beam line, 
allowed for maximum neutron polarization and transmis- 
sion efficiency. A vertical spin guiding magnetic field was 
maintained along the beam line by a combination of per- 
manent magnets and iron plates. The magnitude of this 
field was approximately constant along the polarizer and 
over a large part of the condenser, decreasing only at 
the spin-flippers position to create field gradients, nec- 
essary for spin-flipper operation. Two adiabatic radio 
frequency spin-flippers, mounted around the last section 
of the condenser, were used to reverse the orientation of 
the neutron beam polarization at regular time intervals, 
typically every 16 s. 

The experimental area was shielded from the back- 
ground produced in the neutron guide and polarizer by 
a 0.5 m thick concrete wall with inserted boron-lead col- 
limator (Fig. [3]). The 1.3 m long, multi-slit collimator 
defined the beam cross section to 40x150 mm^ at the 
entrance of the Mott polarimeter. In order to minimize 
neutron scattering and capture, the entire beam line, 
from the collimator to the beam dump, was enclosed in 
a chamber lined with ^LiF polymer and filled with pure 
helium at atmospheric pressure. 

The total flux of the collimated beam was typically 
about lO^'^ neutrons/s. The beam divergence was 0.8° 
in the horizontal and 1.5° in the vertical direction. A 
detailed study of the beam polarization at the position 
of the experimental setup was performed with a polar- 
ization analyzer based on the bent supermirrors concept, 
analogous to the one used in the polarizer. The obtained 
results revealed a maximum polarization of 95% in the 
beam center and its strong dependence on the position 
and on the inclination angle with respect to the beam 
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FIG. 3: Schematic top view of the experimental setup. A sam- 
ple projection of an electron V-track event is shown. The 
coordinate system used throughout this paper is indicated. 

axis. This feature hindered a reliable evaluation of the 
average polarization integrated over the whole beam fidu- 
cial volume. The adopted solution was to measure, in 
parallel to the main correlation experiment, and with 
the same beam and detector, also the neutron /3-decay 
asymmetry. As the asymmetry parameter is known with 
a high precision from other experiments, this approach 
allows the extraction of the average beam polarization, 
while automatically accounting for the complicated beam 
phase space and the detector acceptance. The obtained 
results are listed in Table |lTl 

A more detailed description of the design, operation 
and performance of the cold neutron beam line can be 
found in Ref. [391. 



B. Detector setup and performance 

The Mott polarimeter consists of two identical mod- 
ules, arranged symmetrically on both sides of the neutron 
beam (Fig. [3]). The whole structure was mounted inside 
a large volume dipole magnet providing a homogeneous 
vertical holding field of 0.5 mT within the beam fiducial 
volume. Going outwards from the beam, each module 
consists of a multi-wire proportional chamber (MWPC) 
for electron tracking, a removable Mott scatterer (Pb foil) 
and a plastic scintillator hodoscope for electron energy 
measurement. 

The main requirements which shaped the design of the 
MWPC were the minimal energy loss and multiple scat- 
tering of low energy electrons and the possibly small cross 
section for conversion of gamma quanta into electrons, 
which would pose a dangerous background source. To 
fulfil these conditions a unique combination of special 
features was implemented: 

• The readout of anodes and cathodes allowed reduc- 
ing the total thickness of the MWPC by a factor of 
two. 

• Very light gas mixture based on helium, isobutane 
and methylal (90/5/5), which nevertheless assured 
stable working condition at anode voltages of about 
1800 kV. 
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• Thin, 25 /im Ni/Cr (80/20) wires at 5 mm and 2.5 
mm pitch for anodes and cathodes, respectively, 

• Very thin entrance and exit windows made of 2.5 
/^m aluminized Mylar foil. 

Each chamber contained five planes of anodes (horizontal 
wires) and five planes of measuring cathodes (vertical 
wires), with active areas of 50x50 cm^. The distance 
between anodes and cathodes was 4 mm and between the 
consecutive anode planes 16 mm. The average efficiency 
of a single plane was about 98% and 97% for anodes and 
cathodes, respectively. 

The time measurement of individual wire hits with re- 
spect to the reference signal from the scintillator depends 
on the spatial density distribution of the primary ioniza- 
tion. Unlike in drift chambers, with electrically separated 
drift cells and a sense wire in their center, the adopted 
MWPC geometry does not allow the use of the drift time 
of primary ionization electrons to improve the position 
resolution of a single plane. However, the time informa- 
tion was used to improve the reconstruction of the cluster 
centroid for the cases in which more than two neighboring 
wires have responded. It was also used to check whether 
large clusters are not formed by two overlapping smaller 
clusters. This allowed improving the position resolution 
and the double track resolution of a single plane, and 
was of special importance for cathodes, with an aver- 
age cluster size of about 1.9 hits. The position resolu- 
tion of a single plane, obtained from the distribution of 
reconstruction residua, was about 1.2 mm and 1.7 mm 
r.m.s. for anodes and cathodes, respectively. Better dou- 
ble track and position resolution of anodes lead to more 
precise and more efficient reconstruction of trajectories 
in the vertical coordinate, relevant for the R correlation 
coefficient measurement. 

The scintillator hodoscopes were optimized for the de- 
tection of electrons with energies up to about 1.8 MeV. 
This allowed the distinction of electrons originating from 
neutron decay from more energetic background electrons 
and played a crucial role in the background subtraction 
procedure. Each hodoscope consisted of six 1 cm thick, 
10 cm wide and 63 cm long plastic scintillator slabs. Two 
XP3330 photomultipliers were coupled optically to both 
ends of the scintillator via short light guides of optimized 
shape. This solution allowed for reconstruction of the to- 
tal electron energy with 33 keV resolution at 500 keV 
and 48 keV at 976 keV (see Fig. [S]). The asymmetry of 
the light signal collected at both ends of the scintilla- 
tor slab allowed the determination of the vertical (y) hit 
position with a resolution of about 6 cm, while the seg- 
mentation of the hodoscope in the horizontal direction 
provided a crude estimate of the z-coordinate. Match- 
ing the information from the precise track reconstruction 
in the MWPC with that from the scintillator hodoscope 
considerably reduced the background and random coin- 
cidences. 

Fast pulses from the hodoscope were also used in a 
trigger logic and provided the time reference signal for 



the MWPC wire readout. 



C. Mott scatterer 

Scattering in the field of a spin-less nucleus of electrons 
polarized in the direction perpendicular to the scattering 
plane reveals a left-right asymmetry due to the spin-orbit 
term present in the interaction potential. The purely 
electromagnetic nature of this process (effects of the weak 
neutral currents can be safely neglected in this case) guar- 
anties the exclusive sensitivity to the transverse polariza- 
tion of the incoming electron. The resulting asymmetry 
is proportional to a product of the transverse polariza- 
tion component of the incident electron beam and the 
target-specific analyzing power of the scatterer. The ex- 
periment presented here exploits particularly favorable 
conditions existing for electron scattering on high Z nu- 
clei (lead) at large backward angle, where the Sherman 
function (analyzing power of a single nucleus (40 [) reaches 
its highest value, and the Mott scattering cross section is 
still appreciable (Fig. |4]). 
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FIG. 4: Angular distributions of the Mott scattering cross 
section (left panel) and the Sherman function (right panel) 
for electrons on natural lead. Shaded areas show parameter 
regions relevant to the present experiment. 

For a real scatterer one has to take into account the 
inelastic multiple scattering with atomic electrons and, 
although being much less likely, plural Mott scattering 
at moderate angles which eventually can mimic a sin- 
gle Mott-backscattering event. These effects can signifi- 
cantly deteriorate the initial analyzing power and affect 
the data by their substantial dependence on the thickness 
of the scatterer (c.f. Fig. [5]) and on the incidence angle 
with respect to the foil surface. 

In order to obtain the effective analyzing power of the 
scatterer, Monte-Carlo simulations were performed us- 
ing the Geant 4 simulation framework (41i] and follow- 
ing guidelines presented in Refs. [12, S^. This approach 
takes advantage from accurate theoretical calculations of 
the Sherman functions, properly accounting for atomic 
structure and nuclear size effects, as well as for effects 
due to electrons interaction with the medium. The accu- 
racy of these calculations has been estimated to be better 
than 2% and was verified by comparison with two exper- 
imental data sets: at low (120 keV jH]) and high (14 
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FIG. 5: Simulated effective analyzing power of lead for 
400 keV electrons backscattered at an angle of 140° as a func- 
tion of the lead layer thickness. Electrons incident perpendic- 
ularly to the foil. 



MeV [45|) electron energies. 

In the early phase of this experiment (2003 - 2004) a 
1 /im thick lead layer evaporated on a 2.5 /im thick My- 
lar foil was used as the Mott-scatterer. The bulk of the 
data, however, was collected in 2006 and 2007 with about 
twice as large surface density of the lead scatterer. Even 
this foil was almost transparent to the incident electrons 
from neutron decay, so that more than 99% of them pen- 
etrated through the foil without sizable interaction with 
lead nuclei. 

Due to the Mott-foil manufacturing process, the thick- 
ness of the lead layer was not perfectly uniform. Also the 
illumination of the foil by electrons at the experiment po- 
sition was not uniform. This was included in the system- 
atic uncertainty in the analysis presented in Ref. |28|. In 
order to decrease this uncertainty and to enhance the re- 
liability of the obtained results, precise scans of the lead 
surface density distribution of both 2 /im Mott scatter- 
ers were performed using photon-induced characteristic 
radiation j46| . 

The resulting maps of the lead layer surface density 
measured with an absolute accuracy of about 55 /ig/cm^ 
(Fig.ini), together with the Monte-Carlo simulated multi- 
dimensional effective analyzing power data and the dis- 
tributions of the reconstructed electron vertices on the 
scatterer were used to obtain the final average effective 
analyzing power values. 

The new values are by about 7% smaller than those 
used in Ref. (28j . with relative uncertainty reduced from 
9% to about 3%. 



D. Detector electronics and hardware trigger 

The entire electronics coupled directly to the detectors 
(sense wires of MWPC and photomultipliers) responsible 
for signal amplification, discrimination and derivation of 
the most important early stage of the hardware trigger, 
was designed and built specially for this experiment. The 
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FIG. 6: (Color online) Surface density map of one Mott-target 
with its projection on both y and z axes. The deep minimum 
around y=150 mm is due to imperfect matching of individual 
sheets used for the Mott-target fabrication. Smooth changes 
which are visible in j/-projection are the relics of the evapora- 
tion process used for manufacturing of a single sheet. 



details of the implementation are described in Ref. [47|, 
here only the main concept of the trigger is presented. 

The hardware trigger was built such as to collect vir- 
tually all events belonging to each of two classes: 

• VTl-2 and VT2-1 in Fig. [7]- Mott scattered elec- 
trons with two track segments on one side and one 
segment accompanied by a scintillator hit on the 
opposite side, further referred to as "V-track" (an 
example is shown in Fig. [3]), used for the determi- 
nation of the electron transverse polarization, 

• SI and S2 in Fig. [7] — "single-track" events with 
only one reconstructed track segment on the hit 
scintillator side, used for precise evaluation of the 
average beam polarization. 

In order to enhance the selectivity of the trigger, two 
plane multiplicity signals have been constructed sepa- 
rately for anodes (Y) and cathodes {X) of each detector 
side. High plane multiplicity {XiH, YiH, where i indi- 
cates the side of the detector), was relevant for the de- 
tector side which reconstructed two track segments in the 
MWPC, and low multiplicity {Xn,, Yn,) was required for 
the detector side with only one segment. Taking advan- 
tage of the small drift cell size and using fast OR circuits 
implemented on the discriminator boards of each MWPC 
plane, both signals were generated as early as 80 ns after 
the fastest scintillator pulse. This concept allowed the 
data acquisition to collect V-track events without losses, 
admitting only the most promising candidates for single 
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FIG. 7: (Color online) Simplified diagram of the hardware 
trigger. Shown is the logic of V-track (VTl-2, VT2-1) 
and single track events (SI, S2) constructed from scintilla- 
tor signals (Scl/2) and chamber multiplicities low and high 

1/2 L/H, Yi/2L/h)- 



track events, with the possibility of their further reduc- 
tion by the prescaler module (Fig. [7]). 

Typically the low (high) plane multiplicity signal was 
set when more than two (three) planes of the same elec- 
trode type registered at least one hit wire. With this 
setting the observed rates were about 800 and 8000 Hz 
for V-track and for single track events, respectively. In 
order to keep the dead time at the acceptable level below 
10%, the single track event rate was prescaled by a factor 
of two. 



E. Data acquisition 

For each trigger the information including pulse 
heights and time measurements with respect to the 
fastest hodoscope signal for all hodoscope hits, and 
time measurements for all MWPC wire hits were dig- 
itized in FERA compatible ADC (LeCroy 4300B) and 
TDC (LeCroy 3377) modules. In order to enhance the 
data throughput of the standard FERA bus, a cus- 
tom CAMAC module has been applied (the FERA Tag- 
ger/Extender). It allowed the separation of readout elec- 
tronics into two logical FERA subsystems with separate 
gates [i^. The data from each subsystem were trans- 
ferred in parallel to two pairs of VME hosted memory 
modules working in flip-flop mode. The memories were 
read out via a VME data-bus controlled by a RI02 pro- 
cessor running MBS data acquisition software |^ in- 
stalled on the real-time operating system LynxOS. An 
important role of the data acquisition program was the 
generation of periodic interrupts, typically every second, 
used for read-out of monitoring scalers and for setting 



the spin flippers controlling the beam polarization. This 
software was also responsible for final logging of the data 
on the external mass storage and for sending a fraction 
of the data to the back-end computer for monitoring pur- 
poses. 

The average data flow rate was about 1 MB/s, which 
amounted to about 15% of the maximum achieved data 
throughput of the system. 



IV. 



DATA ANALYSIS 



The final result of this experiment is dominated by the 
last and the longest data collection period. The numbers 
quoted in this chapter apply to this period, however, all 
available data were consistently analyzed in an analogous 
way, and the obtained results are presented and used in 
the calculation of the final average. 



A. Data reduction 

One of the most important features of the data analysis 
is its hierarchical structure. The raw data collected dur- 
ing the experiments, coded in a compact format specific 
to the individual electronic modules, were first converted 
to "physical" format (trajectory segments, deposited en- 
ergy in hit scintillators) with subsequent verification of 
the on-line trigger conditions. All parameters specify- 
ing these conditions were set in such a way that "good" 
events must not be removed regardless of their origin 
{i.e. from- or off-the-beam) . This allowed the prese- 
lection of interesting event classes thereby reducing the 
amount of data to be processed at the next stage of the 
analysis by a factor of about twelve in the case of back- 
scattered events, and by about 50% in the case of single 
track events. 

The second step prepared the final selection of events 
using tighter conditions, thus further reducing the 
amount of data but still allowing for some freedom in set- 
ting the most crucial parameters (listed in Section HV Hp . 



B. Calibration of scintillator hodoscopes 

The energy and position calibration of the hodoscopes 
was performed typically once a week using conversion 
electrons from a ^°^Bi source. A movable support driv- 
ing the source in y and z directions within the symmetry 
plane of the detector {x — 0) between the chambers was 
used to provide uniform illumination of the entire detec- 
tor. 

The reconstructed electron trajectories allowed the 
identification of the hit position along the individual scin- 
tillator and the correction for path-length dependent elec- 
tron energy losses. As a consequence, it was possible to 
calibrate separately relatively short sectors of a scintil- 
lator and to obtain energy calibration specific for the 
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position (y) at which the energy was deposited in the 
scintillator. The reconstructed deposited energy E was 
assumed to be a linear function of the scintillator re- 
sponse function (g): 



E = a,{y) g{Ed,y) + bi{y), i = l 



12. 



(12) 



The function g was defined as the geometrical mean of 
the pulse heights recorded by the "up" and "down" pho- 
tomultipliers (c„, Cd)- 

g = \/Cu Cd- (13) 

With this definition, and assuming uniform light atten- 
uation along the scintillator, the (^-function should be 
proportional to the deposited energy (Ed) and should 
not depend on y. Small deviations from this assump- 
tion were compensated by the position sensitivity of the 
calibration coefficients ai and hi (Fig. [5]). 
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FIG. 9; (Color online) Measured ^"^Bi electron energy spec- 
trum (circles) together with its decomposition into the four 
most important electron conversion lines (482, 558, 976 and 
1052 keV) and a smooth background. The fitted widths cor- 
respond to the energy resolution of the scintillator. 
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FIG. 8: (Color online) Energy calibration for one of the scin- 
tillators (see Eq. (fT2)) '). The residual position dependence of 
calibration coefficients a and h are due to non-uniform light 
attenuation along the scintillator bar. 



The energy resolution can be deduced from the re- 
constructed ^"''Bi electron energy spectrum (Fig. [9]). 
It roughly follows the primary photon statistics and 
amounts to 33 keV around 500 keV. The associated sys- 
tematic uncertainty has been estimated to 5 keV. In order 
to decrease the energy spread caused by different energy 
losses on the way between the electron creation point 
(assumed to be in the symmetry plane of the detector) 
and the scintillator hodoscope, only events with similar 
path-length have been selected for the calibration. 



Using once again the uniform light attenuation as- 
sumption one can show that the signal asymmetry de- 
fined as: 

Eat 



r = In {cu/cd) 



(14) 



should be directly proportional to y and should not de- 
pend on Ed- Figure [TU] presents the y-component of 
hit positions on the hodoscope reconstructed from the 
MWPC information and plotted against the r asymme- 
try. The correlation between both observables is obvious. 
The large width of the band refiects the modest position 
resolution of the scintillator. The average r at a given po- 
sition y was used to obtain the position calibration of an 
individual scintillator y{r). It should be noted that this 
calibration was not perfectly linear, what again indicates 
a small departure from the uniform light attenuation as- 
sumption. 



C. Data selection and event reconstruction 

The first step of the data selection procedure is a fine 
tuning of the coincidence time windows in which ho- 
doscope and MWPC hits are accepted in order to re- 
duce event contamination with accidental coincidences. 
The resulting width of the time window used for the 
MWPC hits (180 ns) accounts for the maximum possi- 
ble drift time of primary ionization electrons in the drift 
cell and for the time walk due to variations of the signal 
rise time. As the time information from scintillators was 
affected only by the time-walk, the length of the corre- 
sponding gate for the hodoscope signals could have been 
much shorter and amounted to 50 ns. 

For a valid hodoscope hit the coincidence between the 
photomultipliers attached to both ends of one scintillator 
was required. In general only one such hit was allowed 
in both hodoscopes. Exception was made for the cases 
when two neighboring scintillator slabs responded. This 
allowed the selection of electrons reaching the hodoscope 
at the scintillators edges which deposited their energy 
partially in two neighboring scintillators. Those events 
played an important role in the determination of the po- 
sition resolution of track reconstruction at the hodoscope 
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FIG. 10: (Color online) A sample position calibration yir) of 
one scintillator bar (solid line) superimposed on the experi- 
mental events distribution. An event is represented by the r 
asymmetry (Eq. I14p and hit position along the scintillator re- 
constructed from the MWPC information. The non-linearity 
due to nonuniform light attenuation along the scintillator is 
clearly visible. Lat corresponds to the average light attenua- 
tion length, Eq. (fTi)) . 



position, based on the position information from the an- 
odes. 

As already mentioned in section fill Bl the information 
from the MWPC was reduced to the cluster centroids. 
The average cluster size was about 1.4 and 1.9 wires for 
anodes and cathodes, respectively. Clusters consisting of 
more than 4 consecutive wires were investigated in order 
to check whether they were formed by two overlapping 
clusters. The splitting condition was derived from the 
analysis of the time information of all hit wires belonging 
to this cluster and relied on the presence of two signifi- 
cant minima in the corresponding time distribution. As 
double clusters are naturally present in vertex topology 
(with their relative distance decreasing for wire planes 
closer to the Mott target), this filter was particularly 
important to increase the detection efficiency of V-track 
events. Clusters larger than 16 wires have been rejected 
as a possible electronics noise (16 channels were grouped 
in one preamplifier-discriminator card). 

Straight lines were then fitted to the obtained hit pat- 
terns separately in anodes and in cathodes using a com- 
binatoric algorithm and the minimum ^ criterion. To 
be accepted a track projection had to be detected in at 
least three wire planes. An event was rejected if at least 
one complete track segment (seen in both projections) 
was registered in excess of the expected number of seg- 
ments, that is exactly one segment at the hit hodoscope 
side and, in the case of V-track events, additionally two 
segments at the opposite detector side. 



In order to reduce background consisting of electrons 
produced in the solid parts of the detector, only track 
segments whose prolongations were contained within the 
active area of the opposite MWPC were accepted. This 
allowed using the two innermost planes of the opposite 
chamber as a veto detector and, as a consequence, to con- 
fine significantly the volume of possible electron origins. 

The reconstructed track segments were confronted 
with a set of conditions checking the consistency of the 
reconstructed event. The extension of the segment recon- 
structed in the MWPC at the hit hodoscope side should 
point, within a given tolerance, to the hit scintillator slab 
and should match the y-position reconstructed from the 
" up-down" asymmetry of the corresponding pulse height 
signals. In the case of V-track events two lines in each 
projection were expected at the detector side which reg- 
istered the Mott scattering vertex. The x-coordinate 
of the vertices reconstructed in both projections should 
match within a tolerance given by the MWPC angular 
resolution. Similarly, matching in both projections was 
required between the segment reconstructed at the hit 
scintillator side and one of the two segments at the op- 
posite side. No scintillator hit in the vicinity of the re- 
constructed vertex was allowed for a valid Mott backscat- 
tering event. 

Additional event classes, absent in the previous publi- 
cation (Ref. [2^) and included in the analysis presented 
in this paper, consist of cases for which the Mott scat- 
tering vertex was detected only in one projection {y-z or 
x-z). Particularly interesting are events scattered close 
to or in the vertical plane. For these the projections of 
both segments in the x-z plane overlap and are recon- 
structed as only one line. The signature of such event is 
weaker: since one vertex is missing one important match- 
ing condition (vertices z-coordinates) drops out. Some 
compensation of this relaxation was achieved by increas- 
ing from three to four the threshold for the minimum 
plane multiplicity of the accepted line fitted in this pro- 
jection. This is justified since each wire along such dou- 
ble track collects twice as much charge as in the normal 
case. It should be noted that the x-z projection was mea- 
sured by cathodes i.e. electrodes with significantly worse 
double track resolution than anodes. This increased the 
number of such events, further referred to as "vertical 
single vertex events" . This event class exhibits maximal 
and exclusive sensitivity to the R correlation, on the av- 
erage two times higher than the primary double-vertex 
event class (see discussion of geometrical form-factors in 
Sect. irVGl) . 

The analogous event class, with one vertex recon- 
structed by cathode planes accompanied by only one 
track segment reconstructed in anodes, is much less nu- 
merous than the vertical single- vertex event class. This 
is due the lower efficiency, position resolution and double 
track resolution of cathodes. Those events are sensitive 
almost exclusively to the N correlation coefficient. 

In the last step of the event reconstruction process each 
event was assigned to one of several event classes deter- 
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mining its role in the further analysis. The most impor- 
tant categories were related to: 

• external conditions: state of beam polarization and 
presence of the Mott target, 

• event geometry: from- and off- the beam, from- and 
off- the Mott-foil. 

The total numbers of reconstructed events are listed in 
Table m separately for all data collection periods. 



TABLE I: Average lead surface density d of the used Mott 
target and total numbers of reconstructed single track events 
S, double vertex events VV and events with one vertex in 
anodes Va and cathodes Vc for all data collection periods. 



by the electron in the magnetic field. Neglecting energy 
losses it can be approximated as: 



Run 


d {fig/ cm ) 


Sx 10^ 


VVxW'^ 


VaxW^ 


VcxW^ 


2003 


1.13 ±0.16 


12000 


19 






2004 


1.13 ±0.16 


43000 


74 






2006 


2.46 ±0.05 


28000 


312 


106 


28 


2007 


2.46 ±0.05 


334000 


1750 


711 


248 


Total 




417000 


2152 


817 


276 



D. Effects of magnetic field 

All matching conditions discussed above as well as the 
event geometry were affected by the magnetic field. 

In order to attain the required uniformity of the spin 
holding field in the beam fiducial volume, the entire de- 
tector was immersed in a constant, large volume mag- 
netic field, produced by a magnet consisting of two soft 
iron plates and eight iron core coils (Fig|3|). Special care 
has been taken in order to shield all photomultipliers 
against the influence of this field. This has been achieved 
with a double layer of mumetal shielding around each 
photomultiplier . 

The effect of the magnetic field on the detected elec- 
trons is twofold. Their spins precess at the Larmor fre- 
quency and their trajectories are bent so that their x-z 
projection becomes an arc. For the measurement of the 
transverse polarization of electrons those effects are po- 
tentially dangerous, since the polarization is produced 
and analyzed at distant locations (neutron decay and 
Mott-scattering) . However, since the electron g- factor 
is almost equal to 2, the spin precession almost exactly 
follows the momentum rotation. The maximum remnant 
effect, due to the "g — 2" factor, is well below one arc 
minute in this experiment and is therefore irrelevant for 
the achieved accuracy. 

There are, however, other consequences of the bend- 
ing of the electron trajectories. Part of them can be ac- 
counted for while others are discussed in order to demon- 
strate why they do not influence the final result. 

The deviation p of the electron trajectory from the 
straight line depends strongly on the distance d traveled 



dtan 



"1 . / 




— arcsm 









(15) 



where R is the curvature radius which, for the lowest 
electron energy detected in this experiment, amounts to 
about 3 m. From this it follows that in the worst case 
the deviation of the real trajectory from a straight line 
within one MWPC is below 0.5 mm. Considering the 
average size of clusters in cathodes (SS14 mm), such a 
small correction can be neglected and one may safely 
use a linear fit to the data. This is not the case if one 
considers prolongations of the obtained lines to the Mott 
target, to the scintillator hodoscope or to the opposite 
MWPC. Then the effects can be substantial, but knowing 
the electron energy and the magnetic field strength they 
can be accounted for. 

A strict correction would require taking into account 
continuous energy losses along the electron path. How- 
ever, in view of the much larger effect of electron multiple 
scattering, this correction has been simplified by using 
the average electron energy along the considered path 
segment. Subsequently the line fitted in the MWPC was 
treated as a tangent to the circular trajectory with a ra- 
dius corresponding to this average energy. Trajectories 
obtained in this way were then used to calculate extrap- 
olated electron positions and incidence angles at beam, 
hodoscope, Mott target, etc. 

Surprisingly, even a weak magnetic field can have a 
significant influence on the efficiency of V-track recon- 
struction. In order to understand this effect two kinds 
of V-tracks must be introduced. In the following they 
will be referred to as convex and concave (Fig. [TT]) . For 

Mott target ® ^ 




FIG. 11: Convex and concave vertices. 5 corresponds to the 
displacement of the electron trajectory due to the magnetic 
field, measured along the z-axis at each wire plane. 

convex V-tracks the distance between both trajectories 
measured within each wire plane (along the z direction) 
is always larger than for its concave analog. This differ- 
ence reaches its maximum around the outermost planes, 
closest to the Mott target, at the place where the sepa- 
ration of clusters belonging to both arms of the vertex 
reaches its minimum. This causes a difference in the 
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probability for two clusters in the same plane to over- 
lap and, as a consequence, increases the detection effi- 
ciency of convex V-tracks and decreases this efficiency 
for concave V-tracks. Of course the same effect, however 
with decreasing significance, occurs also in other MWPC 
planes. 

Another effect may be described as focusing (for con- 
vex) or defocusing (for concave) of the long arm of a 
V-track on the active area of the opposite detector. It 
acts coherently with the previous one, further increasing 
(decreasing) the detection efficiency for convex (concave) 
V-tracks. 

In order to reduce the impact of all effects induced 
by the guiding magnetic field, its magnitude has been 
reduced from 1 mT, used between 2003 and 2006, to 
0.45 mT in the 2007 data taking period. Nevertheless 
the effects persist and can readily be observed (see e.g. 

Figs, na HI. 

E. Background correction 



of the tracks at the opposite detector side [32|. These 
conditions have to account for the angular resolution of 
the MWPC, the beam density distribution and its diver- 
gence. The better signal to background ratio obtained 
for single track events can be attributed to their much 
larger number what allowed for a much tighter setting of 
all geometrical cuts. 

The validity of this background subtraction method 
was verified by comparing background-corrected energy 
spectra with the simulated /3-decay spectra in which en- 
ergy losses and detector resolution were taken into ac- 
count. 
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Two kinds of background have been taken into account 
and corrected for. The first, further referred to as the 
"off-beam" background is present in single track and in 
V-track event classes, while the second, "foil-out" back- 
ground, applies only to the V-track events. 

For the "off-beam" background the number of electrons 
not originating from the free neutron decay, was deter- 
mined by comparing energy spectra of two event classes: 
(i) events for which the reconstructed electron trajectory 
crossed the neutron beam volume ( "from beam" ) and (ii) 
events for which the electron origin was outside the neu- 
tron beam ("off beam"). The procedure relies on the 
assumption that the spectral shape of the background is 
the same for both event classes, while the characteristic 
neutron /3-decay spectrum with end-point energy of 782 
keV is present only in the "from beam" class (Fig. [T^ . 
This allows scaling the "off beam" background distribu- 
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FIG. 12: (Color online) Energy distributions of signal (green) 
and background contributions (red) for single track events 
(left) and double vertex event class (right). Siiaded areas in- 
dicate the "signal" energy range used in the final calculations. 

tion such that it matches the high energy part of the 
"from beam" energy spectrum. For this assumption to 
hold, the "off beam" range has to be carefully chosen for 
both: the inclination angle and the extrapolated origin 



FIG. 13: (Color online) Background-corrected experimental 
energy distributions (shaded areas) of (a) single-track and (b) 
double-vertex V-track events compared with simulations. 

Such a comparison is shown in Fig. [T^] for single tracks 
and for the Mott scattering events. In the latter case 
the modification of the /? spectrum induced by the en- 
ergy and angular dependence of the Mott-scattering cross 
section is clearly visible. Electronic thresholds are not in- 
cluded in the simulation — this is why the measured and 
simulated distributions do not match at the low energy 
side. In contrast, the matching at the high energy side 
is nearly perfect for single track events (Fig. [T4|) . Due 
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FIG. 14: Detailed view of the energy distributions of the 
"from beam" and "off beam" single track events in the vicinity 
of the electron end-point energy for a defined range of elec- 
tron emission angles. The agreement between the measured 
(solid line) and simulated signal distribution (full circles) is 
near to perfect. 

to not well defined energy losses (in particular the deter- 
mination of the depth at which the Mott scattering took 
place within the lead scatterer is far beyond the accuracy 
of the electron tracking) the analogous comparison, at a 
similar level of accuracy, is not possible in the case of 
V-track events. 
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The lack of two important matching conditions is the 
reason of the worse signal to background ratio observed 
in event classes with the single vertex signature fFig. lTSi) . 
The difference between events with the vertex recon- 
structed either in cathodes or in anodes can be explained 
by the worse double track resolution of cathodes and the 
horizontal geometry of the beam. 
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FIG. 15: (Color online) Energy distribution of signal (green) 
and background contributions (red) for single vertex events 
in anodes and cathodes (left and right, respectively). Shaded 
areas indicate the "signal" energy range used in the final cal- 
culations. 



in the off-the-beam volume with more obtuse vertices. 

The above argument does not apply to events with 
a single vertex in cathodes (Fig. [T71). In this case the 
shape of the background resembles that of single tracks 
but the signal to background ratio is much worse than 
for events for which only one vertex was reconstructed 
in anodes. This can be explained by the beam fiducial 
volume geometry which is well defined in the vertical di- 
rection (±90 mm) but with horizontal limits fixed only 
as a compromise between striving to enhance statistics 
and to minimize the background of electrons originating 
from the MWPC frames. 
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The same background type can be observed in the y- 
projection of the position distributions of the extrapo- 
lated origins for low and high energy electrons. The ori- 
gin of the electron track was taken to be the intersec- 
tion between its trajectory and the symmetry plane of 
the detector (x = 0). This approximation is justified by 
the narrow beam size along the x-coordinate (± 2 cm) 
defined by the ®Li beam collimator. The clear profile 
of the neutron beam can be recognized in the position 
distributions of low energy electrons (E < 750 keV) in 
contrast to the distributions of events with higher energy 
(Fig. [TB|) . The significant difference in the background 
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FIG. 16: (Color online) Beam profile seen in the position dis- 
tribution of the extrapolated electron origins for single track 
(left) and V-track events (right), for low and high electron 
energies. Shaded areas indicate the position range accepted 
as "from beam" region. 

distribution for double vertex events as compared to the 
case of single track events can be explained by the tighter 
setting of the geometrical cuts and by the dependence of 
the V-track detection efficiency on the Mott scattering 
angle. As already discussed, for more acute vertices the 
probability that clusters overlap increases, what as a con- 
sequence enhances the contribution of events originating 



FIG. 17: (Color online) Beam profile seen in the position dis- 
tribution of the extrapolated electron origins for events with 
a single vertex in anodes (left) and cathodes (right), for low 
and high electron energies. Shaded areas indicate the position 
range accepted as "from beam" region. 

In the case of V-track events, beside the background 
discussed above, events for which the backscattering took 
place in the surrounding of the Mott target induce an ad- 
ditional source of background. Figure [T8l presents the dis- 
tributions of the reconstructed x-component of the vertex 
positions for the data collected with and without Mott 
target. The "foil-out" distributions have been scaled ap- 
propriately by a factor deduced from the accumulated 
neutron beam for each setting of the Mott scatterer. 

The "foil-in" distributions clearly peak at the Mott 
target position. The broad maxima observed in "foil- 
out" distributions can be explained by backscattering on 
the MWPC material (mainly on the Ni-Cr wires and the 
aluminized exit window) and on the wrapping of scintil- 
lator hodoscopes. The very good signal to background 
ratio observed for the double vertex events (« 23) de- 
creases for single vertex events in anodes (w 16) to reach 
its minimum for events with a single vertex in cathodes 
(« 9). The sharp structures seen in the distributions of 
one-vertex events, for |a:;| < 200 are due to an artifact of 
the reconstruction procedure caused by the assignment 
of single wire clusters to the discrete wire positions. 



F. Beta decay asymmetry 

To extract the beam polarization P averaged over the 
beam fiducial volume the following asymmetries were an- 
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FIG. 18: (Color online) "Foil-out" background contributions 
(red) to the vertex a;-coordinate position distributions of V- 
track events, for the left and right detector side (left and 
right panels, respectively), separately for double- vertex events 
(top panels) and events with only one vertex (reconstructed 
in anodes - middle, and cathodes - bottom panels). The 
arrows show the Mott foil position. Shaded areas indicate the 
position range accepted as "from-foil" region. 



alyzed: 



f (/3,7) 



iV+(A7)-^-(/?,7) 
iV+(/3,7)+iV-(/3,7)' 



(16) 



where are experimental, background corrected num- 
bers of counts of single tracks, sorted in 4 bins of the elec- 
tron velocity normalized to the speed of light (3 = v/c, 
and 11 bins of the electron emission angle 7 with respect 
to the neutron polarization direction. Considering only 
the relevant terms in Eq. ([2]), can be written as: 



iV± (/3, 7) = iVo [1 + ry* A P ■ /3 JF(/3, 7)] , (17) 

where the sign in superscripts reflects the beam polariza- 
tion direction and /3J-" is a kinematical factor correspond- 
ing to the average z-component of the electron velocity 
in a given bin of /3 and 7: 



/3^(/3,7) = 



•J) 



/3,7 



(18) 



with J a unit vector in the direction of the neutron po- 
larization. The factors rj^ account for the spin flipper 
efficiency. When the spin flipper is switched off one has 
r]- = l (original polarization is fully maintained) , whereas 
in the opposite case 77+ = — 1 -f 77. The value w 0.0114, 
indicates that only a very small fraction of all neutrons 
had not reversed their spin after passage through the spin 
flipper generated RF field. The magnitude of this effect 
has been investigated in a dedicated experiment ^ . The 
factors account for the influence of the spin flipper op- 
eration on the detection efficiency of single track events. 
The radio frequency associated with the spin flipper op- 
eration, propagating via electrical grounding and in the 



air via the beam line volume, increases slightly the noise 
level observed on the wires of MWPCs and in the ho- 
doscopes. This increases the dead time and decreases 
the reconstruction efficiency. Similarly, as in the previ- 
ous case = 1 and e"*" = 1 — e is close to, but less than 
unity. The actual value of e (sa 0.004) can be calculated 
from the total numbers of single track events accumu- 
lated in each beam polarization state and corrected for 
the corresponding beam intensity. It is visible as a small 
negative offset of all experimental points in Fig. [111 With 
the above definitions, applying first order Taylor expan- 
sion in e, 7] and second order in the term P A f3 T{l3,j), 
Eq. ([HI) reads: 



£:(/?, 7) = (l-77/2)PA/3^(/3,7)-e/2. (19) 

Figure [THj shows the obtained £ as a function of /3 F for 
different electron energy ranges. Taking A as a constant 
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FIG. 19: (Color online) Asymmetry E (Eq.[l9l) of signal (full 
symbols) and background (open symbols) as a function of 
/3 F for different electron energy ranges for events registered 
on one detector side. The fit of Eq. (|19|l to the data allows 
to extract the neutron beam polarization. A tiny negative 
shift of all E values is due to the spin flipper influence on the 
detection efficiency e^. 

known with a very good precision (^A = —0.1173± 0.0013 
50]), the average neutron polarization can be obtained 
from a one-parameter fit of Eq. (|19p to the experimental 
data. 

It should be noted that in the background correction 
procedure special care has been taken to ensure that the 
background counts do not depend on the beam polar- 
ization direction. The E asymmetries for events which 
do not originate from the beam volume (Fig. 1191 open 
symbols) as well as of the high energy events (above the 
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electron end point energy) do not depend on the angle 7 
and are consistent with zero polarization of their sources. 

The average neutron polarization values for the four 
data taking periods are listed in Table |lll The low po- 
larization for the 2004 data set has been traced to a bug 
in the guiding field found post factum and verified in a 
dedicated experiment. 

The electron emission asymmetry should also be ob- 
served in Mott scattered event classes. Due to the much 
lower statistics of those events, the extracted polarization 
Py is much less precise. However, a similar analysis has 
been performed and its results are in satisfactory agree- 
ment with the single track data (Table |lT| . 

TABLE II: Summary of neutron beam polarization analysis. 
Polarizations deduced from double and single vertex event 
classes {Pvv and Pv) are also shown. 



Run 


P X 10^ 


Pvv X 10^ 


Pv X 10^ 


2003 


80.3±1.3±1.6 


71.8±9.4±1.6 




2004 


44.2±0.4±1.5 


48.7±8.3±1.5 




2006 


80.0±1.0±1.5 


82.9±3.9±1.5 


74.1±9.5±1.5 


2007 


77.4±0.2±0.7 


78.7±1.7±1.2 


79.9±3.0±1.2 



For the longest data taking period in 2007, an inde- 
pendent analysis of the beam polarization has been per- 
formed. This approach used the same asymmetries as 
defined in Eq. (|16p. however, with different binning. The 
single track events were sorted in only 2 bins in 7 (elec- 
tron emission into lower and upper hemispheres) but, ad- 
ditionally, in 2787 time bins (half an hour long each). 
This allowed to search for a possible time dependence of 
the extracted polarization and daily modulations of dif- 
ferent observables, and also provided a consistency check 
with the previous analysis. The constant value fit to the 




FIG. 20: Time series of polarization extracted from the analy- 
sis of single tracks sorted in 2787 time bins. The longer breaks 
correspond to no-beam periods. 

data fFig. [20]) results in an average polarization of 0.773 
± 0.002 which is in a very good agreement with the result 
of the previous method. Considering only the statistical 
errors, the P^r degree of freedom of this fit amounts 
to 0.94. The discrete Fourier transform of the obtained 
time series is consistent with a white noise distribution. 
One can conclude that no significant time dependence has 
been observed within about 3 months of the 2007 data 



collection period. For more details of this analysis and 
obtained limitations on some Lorenz invariance violating 
parameters see Ref. [sij . 



G. Correlation coefficients R, N 

For the analysis of the transverse electron polarization 
components and the associated correlation coefficients, 
the method of backward Mott scattering was applied. 
Parity and time reversal conservation of the spin-orbit 
force responsible for the spin dependence of this electro- 
magnetic process guarantee its exclusive sensitivity to the 
transverse polarization component perpendicular to the 
scattering plane. Technically, this can be expressed by 
the following substitution: 



(20) 



where S is the effective analyzing power of the Mott- 
scatterer which for a given target element depends on 
electron energy E and scattering angle 9, and fi is a 
unit vector perpendicular to the Mott-scattering plane 

(Fig.m 



Pe X Ps 
|Pe X Ps 



(21) 



where Pe and Ps are incident and scattered electron mo- 
menta, respectively. 




FIG. 21: Definition of angles relevant for the analysis of the 
transverse electron polarization by Mott scattering. J is the 
neutron spin and Pe, Ps are incident and scattered electron 
momenta, respectively. 

Applying this substitution to Eq. ^ , the background- 
corrected experimental numbers of counts of V-track 
events can be expressed as: 



no e 



A/3 J'(a)- 



NSgia) + RpSn{a) 



(22) 



The sign in superscripts reflects the beam polarization 
direction, the meaning of and 77^ is the same as in the 
case of single track events, Eq. (|T7)) . and the kinematic 
factors ^(a), g{a) and 'H{a) represent the average values 
of the quantities J • p, 3 ■ a and J • p x &, respectively 
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(Fig. [22]) . The bar over a term indicates event-by-event 
averaging used in all analyzing methods applied in this 
work. In order to fully exploit the symmetry properties 
of both the physical problem and the experimental setup, 
all those quantities where sorted in 12 bins of a, defined 
as the angle between the electron-scattering and neutron- 
decay planes (Fig. [^T|) . represented respectively by the fi 
and rh unit vectors. 



J X Pe 
|J X Pe 



(23) 
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FIG. 22: (Color online) Average geometrical {P, Q, Ti.) and 
kinematical (/3, S) factors as a function of a for double (open 
symbols) and single vertex events (full black and red symbols 
refer to vertical and horizontal V-tracks, respectively). Dot- 
ted lines are to guide the eye only. Error bars are smaller 
than the symbol size. 

The a angle distribution of the sum n+ + for all 
classes of the analyzed V-tracks is shown in Fig. [221 The 
deep minima around integer multiplicities of 7r/2, in the 
distribution of double vertex events (VV) reflect the rect- 
angular geometry of the MWPC and are mainly shaped 
by limited double track resolution of anodes and cath- 
odes. The significantly lower intensity at negative values 
of a angle (corresponding to concave V-tracks) is due 
to the magnetic field infiuence on the V-track detection 
efficiency, discussed in Section HVDl 

The visible spread of the + distribution for "ver- 
tical" single vertex events [Va) is entirely an effect of the 
energy dependent correction for the magnetic field ap- 
plied to the electron trajectories. Without this correction 
only two discrete values (0 and tt) would be possible (all 
relevant vectors Pe, Ps and J are coplanar in this case). 

Two different approaches have been used to obtain the 
N and R correlation coefficients. The first one, presented 



10= 



10' 



10^ 



10= 



-2 2 

a 



FIG. 23: (Color online) The a distribution (n+ +n~) of all 
analyzed double vertex events ( W) and events with only one 
vertex reconstructed in anodes (Va), and cathodes (Vc). 



in the next section, can be applied to both V-track event 
classes, those with full geometrical information and those 
with only one vertex in the vertical plane. The results of 
this approach have been adopted as the final result of this 
experiment. The second approach allows the extraction 
of N and R coefficients separately from dedicated double 
ratios. However, in the case of the R correlation this 
approach requires the assumption that the experimental 
data are symmetric with respect to the transformation 
a a' = —a. This requirement is drastically violated by 
the influence of the spin holding magnetic field in the case 
of the "vertical" V-track event class. In the case of events 
with full geometrical information the a — > —a symmetry 
is much better fulfilled, allowing for the final application 
of the double ratio method, albeit, as a consistency check 
only. 



1. Correlation coefficients R, N from asymmetry 

To extract the N and R correlation coefficients the 
following set of asymmetries was considered: 



A{a) 



Ti^ (a) — n {a) 
n+ (a) -|- W {a) 



(24) 



Applying Eq. 



and the first order Taylor expansion 



in small quantities (ry, e^, RSjST-L and NSQ) and second 
order Taylor expansion in the largest term, PA(3T, one 
obtains: 



Aia) = Pil-T]/2) A(3J'{a) + 



NS{a)g{a)+RSia)pn{a) -fi,,/2, (25) 

The term PA(iF accounts for the nonuniform illumina- 
tion of the Mott foil due to the /3-decay asymmetry and 
is known precisely from event- by- event averaging. The 
systematic uncertainty of this term is dominated by the 
error of the average beam polarization P. It is interesting 
to note that the functions Q and Ti. follow quite closely 
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FIG. 24: (Color online) Experimental asymmetries A cor- 
rected for the PApj- term as a function of a. The solid line 
illustrates the two-parameter (A^, R) least-square fit to the 
data using experimental form factors SQ and j3S'H. The in- 
dicated errors are of statistical nature. 
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FIG. 25: A'^ correlation coefficients calculated according to 
Eq. (|28|l for double vertex events (full symbols) as a function 
of a and their average. The indicated errors are of statistical 
nature. Open symbols: the same but for unpolarized beam 
(see Sec. lIVG3|) . 



sine and cosine functions, respectively, and are almost 
orthogonal to each other. As a consequence, the covari- 
ance matrix of the two-parameter fit used to obtain the 
A'^ and R correlation coefficients is almost diagonal, with 
the correlation coefficient p{R, N) w0.007. 

A two parameter fit of the experimental asymmetries 
A, corrected for the PA(3T term, to the experimental 
data set of 2007 is shown in Fig. [Ml The extracted values 
for the R and N coefficients are listed in Table [ 



2. Correlation coefficients R, N from double ratios 

From the approximate symmetry of the detector with 
respect to the transformation a — > — a, it follows that /3, 
S and the factors J^, % are almost symmetric while Q 
is an almost antisymmetric function of a (Fig. [22]) . Ap- 
plying these symmetries, the Taylor expansion as in the 
previous section and the definition of the quantities n^, 
Eq. (|22|) . one can see that the double ratio defined as: 



Q{a) = 



(r(a) + l)' 



where 



r{a) 



'nr^{a)n (—a) 
n^{a) ri+(— a) ' 



(26) 



(27) 



allows to extract the N correlation coefficient according 
to: 



l-l{(l-j^/2)PApFy 
(1 - Ti/2)PSg 



(28) 



The advantage of this method is that the effect associ- 
ated with the term PA/ST is suppressed by a factor of 
about 60 as compared to Eq. The ratio Q is also 

insensitive to the spin fiipper related modulation of the 



detection efficiency. Figure [25] shows the values of N ob- 
tained as a function of the angle a with their average 
value. The good agreement between the N values ob- 
tained in both ways (Table IIIII) enhances our confidence 
in the experimental values of the N and R coefficients 
obtained in the previous section. 

An alternative way to extract the R correlation coeffi- 
cient makes use of the analysis of another ratio: 

r'(a)~l 



where: 



r'{a) 



In (a) n {—a) 
71+ (a) n+(— a) 



(30) 



Applying Eq. (|22[) and keeping only terms linear in small 
quantities (e^, , PA/3F, PR/3SQ), one can show that: 



R 



U -il-fj/2)PA/3F-e,./2 



(31) 



{l-T]/2)P/3SH 

In this method one suppresses the term proportional to 
the N correlation. It is, however, sensitive to the "false" 
asymmetry due to the term PA^T and to the spin flip- 
per related modulation of the detection efficiency, and 
therefore has no clear advantage over the method based 
on Eq. (|25p. Moreover, since the distribution of events 
with only one vertex in anodes is not symmetric with 
respect to the transformation a — ^ —a (Fig. [23[) . this ap- 
proach can not be applied to this event class. Figure [^Sj 
shows the obtained values of i? as a function of the angle 
a with their average. The results for the double vertex 
event class are shown in Fig. [26] as a function of a with 
their average, and are also included in Table [Till 

3. Polarization of background and unpolarized beam 

An important consistency check of the analysis of the 
Mott scattered events relied on the determination of the 
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TABLE III: Summary of results obtained in all data collection periods. Statistical and systematic uncertainties follow the 
experimental values. Nsm is the SM value of the A'' coefficient calculated at Ek- Its error comes from the experimental 
uncertainty of the decay asymmetry parameter A [H^l . The difference as compared to the results presented in [2^ is due to 
an improved determination of the effective analyzing powers and including the additional class of single vertex events (V^). 
/n.d.f. of all fits is also presented. 



Run Ev. Class iVsM x 10^ A^xlO^ 






iVxlO^ 


xVn.d.f. 




xVn.d.f. 


Eq. ^ 


Eq. 


Eq. (123 


Eq. 


Eq. ^ 


Eq. dsn 


Eq. 



2003 


VV 


71±1 


89±92±31 


-90±137±38 


1.6 


139±124±27 


1.9 


-55±152±42 


1.4 


2004 


W 


68±1 


74±80±17 


-135±130±30 


1.8 


171±103±15 


2.0 


-58±148±30 


1.5 


2006 


W 


68±1 


94±35±10 


-13±48±10 


1.3 


97±35±10 


0.6 


-36±48±12 


2.1 


2006 


V 


68±1 


44±109±23 


-50±55±21 




53±117±23 








2007 


W 


68±1 


59±13±5 


13±18±6 


1.1 


63±14±5 


1.3 


-5±18±6 


0.7 


2007 


V 


68±1 


51±32±14 


9±20±13 




52±33±5 








Total 






62±12±4 


4±12±5 




67±11±4 









0.1 



QC 



-0.1 









R = -0.( 


)05 + 0.01 


8 












> 


1 





7l/2 

a 



FIG. 26: R correlation coefficient calculated according to 
Eq. pi[) for each a bin for double vertex events. Only statis- 
tical errors are indicated. 



background polarization. As there was no conceivable 
mechanism which could cause such polarization the ex- 
pected value was zero. To check this presumption, the 
asymmetries as defined in Eqs. ([T6l) and (f24| have been 
calculated for events with energy larger than the neutron 
/3-decay end-point energy (Figs. [27] and [29l) and for events 
originating outside the beam fiducial volume (Fig. I19p . 
It turned out that within the statistical accuracy all 
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FIG. 27: Asymmetry £ (/3, 7) for single-track events with en- 
ergy larger than the neutron /3-decay end-point energy, sepa- 
rately for the left and right detector side. 
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FIG. 28: Asymmetry £ (/3, 7) for single-track events and un- 
polarized beam. Both detector sides have been added to- 
gether. 
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FIG. 29: Asymmetry A{a) for double-vertex events with en- 
ergy larger than the neutron /3-decay end-point energy (open 
symbols), and for unpolarized beam (full symbols). 



asymmetries were consistent with zero. This proves that 
the analysis was not biased by, for instance, a spin flipper 
related false asymmetry. 

The same test has been performed for the data taken 
with an unpolarized beam. In this case, the analysis dif- 
fered from the regular analysis of polarized beam in only 
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one detail: each two consecutive spin states correspond- 
ing to one flipper "on" and one flipper "off" periods were 
concatenated into a new one called "spin-up" , while the 
next two were used to obtain a new "spin-down" state. 
This construction assures averaging of the polarization 
over an equal number of neutrons in both spin states and 
hence leads to an unpolarized beam (we neglect here the 
effect due to spin flipper inefficiency which is well below 
the statistical accuracy of this analysis). 

Also in this case the asymmetries were consistent with 
zero (Fig. [25)1 what strengthens the confidence in the ob- 
tained final results. 



H. Systematic errors 

The systematic uncertainties involved in the evalua- 
tion of the R and N coefficients are dominated by effects 
introduced by the background subtraction procedure and 
the choice of specific values of the cuts which determine 
whether an individual event is attributed to the "signal" 
or to the "background" . 

A "signal" event is defined as an electron originating 
from the free neutron decay and backscattered off the 
Mott target. From this definition it follows that the set 
of necessary conditions describing such event must in- 
clude: (i) geometrical limitations to the beam volume, 
(ii) specification of the allowed energy range and (iii) geo- 
metrical limitation to the area of the Mott target applied 
to the reconstructed scattering vertex. As a general rule, 
the symmetry of the detector setup has been preserved 
in the definition of cuts. The single exception from this 
rule was made for the beam limitation along the y coor- 
dinate: the cutting line was inclined by an angle of about 
1° in accordance with the beam divergence. This allowed 
limiting the number of parameters to the following three 
groups (see also Fig. [S]): 

• "from/off beam" definition: zi^^_^, ?/2„,„^, 

• from/off neutron decay: E^^.^, -E'i,„ax: -^^fminJ 

• from Mott foil: Xmin, Xmax, Wmax (sidc length 
of the square indicating the foil area in the y—z 
plane) . 

Those parameters were used to classify an event as be- 
longing to: 



1. signal: 



\y\ < Vlr^a.. \A < (32) 

> X > Xmax, \Y\ < Wmax, \Z\ < Wmax, 



2. off beam background: as above but 

\y\ > y2„„. (33) 



3. high energy background: 

< yi,„ax, kl < ^Ir^a., (34) 

E > En^i^, E < Eh^^^i 

Xmin > X > Xmax, \Y\ < Wmax, \Z\ < Wmax, 

where (a; = 0, y, z) and {X, Y, Z) denote coordinates of the 
electron origin and the Mott scattering vertex, respec- 
tively. Due to the limited accuracy of the reconstructed 
energies and trajectories, for each of those parameters 
there exists a certain range of values which seems to be 
almost equivalent. In the analysis, however, the specific 
values can generate slightly different final results. In or- 
der to estimate this effect, the corresponding ranges of 
acceptance have been identified individually for each pa- 
rameter and the final analysis was repeated varying one 
parameter in its range with all others fixed in the center 
of their ranges. The maximal deviations of the resulting 
R and N coefficients from the central values were taken 
as maximal errors. Figure [501 presents the result of such 
analysis for parameters limiting the geometrical position 
of the Mott scattering vertices. The contributions of all 
parameters to the final uncertainty, as well as their ranges 
used in this analysis, are collected in Table HVl 
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FIG. 30: Influence of the geometrical definition of the ac- 
cepted scattering vertex positions on the reconstructed R and 
A'' correlation coefficients. Vertical lines indicate values ap- 
plied in the analysis. 

It should be noted that the vertical alignment of the 
apparatus with respect to the neutron beam has been ver- 
ified to a precision below 1 mm using the reconstructed 
centroid of the beam profile fFigs.[T6land [T7|) . A similar 
precision of alignment with respect to the magnet setup 
and the beam axis was maintained for all other detector 
components except the Mott scattering target. Here the 
accuracy of positioning in the a;-coordinate was about 
±2 mm. This is still acceptable considering the vertex 
reconstruction accuracy of the MWPCs (Fig. [T8|. 

The next important systematic uncertainty is due to 
the limited accuracy of the determination of the average 
beam polarization. All measured asymmetries used for 
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the evaluation of individual correlation coefficients are 
proportional to the product of this coefficient and the 
beam polarization. As a consequence, the relative error 
of the extracted coefficient must be larger than that of 
the polarization. 

The situation is even more difficult for the R correla- 
tion. With a vertically polarized neutron beam the exis- 
tence of a nonzero value of this correlation would result 
in a difference between the number of electrons backscat- 
tered into the upper and into the lower hemisphere. The 
same effect is generated by the /3-decay asymmetry {A 
correlation), appearing as a term PA/3T in Eqs. ([25)1 
and ((3T|) . What makes the link between the A and R 
correlations so special is the shape of the corresponding 
average form factors T and H (Fig. [22]) . They both ex- 
hibit the same symmetry properties with respect to a, 
so that the effects generated by each of them are almost 
indistinguishable. 

The cleanest, but also unpractical, way to avoid this 
interference would be to confine the electron emission an- 
gle to 90°, in which case .F = 0. With a ffirite accepted 
solid angle one is forced to apply a suitable correction in 
the form of the PA(3J- term. Since the electron momenta 
after emission and after scattering are reconstructed, the 
necessary form-factors {F, G, H) are known to a high pre- 
cision for the entire event sample, and the main impact 
of the applied correction is due to the uncertainty of the 
neutron polarization. To evaluate the magnitude of the 
influence of this term on the final result, the fit with R 
and N as free parameters was repeated with P varied by 
one standard deviation of the total neutron polarization 
uncertainty (Fig. [3T|) . The obtained difference enters the 
budget of the systematic errors and is presented in Table 



neutron polarization direction. A residual systematic ef- 
fect (Table IIV|) was induced by the uncertainty of the 
field measurements and by the observed stability of the 
ambient magnetic field at the experimental position. 

Since the radio-frequency signal of the spin flipper was 
a source of small noise in the readout electronics, tiny 
dead time variations correlated with the spin flipper were 
observed. Their influence on the result was corrected for. 
The residual effect is included in Table ITVl 

TABLE IV: Summary of systematic errors for the 2007 data 
set. The ranges over which individual parameters have been 
varied in the error estimation procedure are shown. 

Source 5N x 10* 5R x lO" 



Total 54 61 



In the final error analysis it has been assumed that 
the sources contributed independently, and so these were 
added quadratically to obtain the final systematic uncer- 
tainty. 
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FIG. 31: Dependence of the reconstructed 7? and A'^ correla- 
tion coefficients on the value of the neutron beam polariza- 
tion. Vertical lines indicate values applied in the analysis. 
The strong correlation between R and P is due to the PpAJ" 
term in Eq. (f25|) . 

Despite the careful design and manufacturing of the 
large volume electromagnet responsible for the spin hold- 
ing magnetic field, the mapping of this field at the ex- 
perimental position showed a small misalignment with 
respect to the vertical direction and nonuniformities in 
the beam fiducial volume. The average effect was ac- 
counted for in the analysis by appropriate rotation of the 



V. RESULTS 

Combining the results from all runs leads to the final 
result (Table mi): 

R = 0.004 ± 0.012,tat ± 0.005«yst, (35) 
N = 0.067 ± O.Oll^tat ± 0.004,yst. (36) 

In Figs. |32] and |33] the new results are presented as 
exclusion plots containing in addition the experimental 
information available to date from nuclear and neutron 
beta decays, as surveyed in Ref. [U]. The upper part 
of Fig. 131] contains plots corresponding to the real and 
imaginary parts of the normalized scalar and tensor cou- 
pling constants S and T, Eqs. (jH]) and (ITUl) . The present 
accuracy of the determination of the N correlation co- 
efficient does not improve the already strong constraints 
on the real part of the couplings (left panel). It is, how- 
ever, consistent with the existing data and, in addition, 
adds confidence to the validity of the extraction of the 
R correlation coefficient. The latter constrains signifi- 
cantly the imaginary part of the scalar couplings, be- 
yond the limits from all previous measurements (right 
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panel). Also in this case the result is consistent with the 
SM (Cs = = Ct = = 0) and with time-reversal 
invariance. 

In the lower part of Fig. [32] the same convention was 
used to illustrate the constraints (existing and resulting 
from the present work) to the helicity projection ampli- 
tudes in the leptoquark exchange model, as defined in 
Ref. [s^l. In this formalism FlLj/lh and ^^LL,ft.Lfl 
correspond to leptoquarks with charge \Q\ = 2/3 and 
IQI = 1/3, respectively. Capital letters correspond to the 
scalar (spin-zero) while lower-case letters describe vector 
(spin-one) leptoquarks exchange amplitudes. Subscripts 
indicate the helicity structure of the underlying interac- 
tion. As in the previous case, only the R correlation 
reveals evident exclusion power and allows a significant 
improvement of the constraints on the imaginary part of 
the vector leptoquark amplitudes {/lr. + hLn)- 
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FIG. 32: Experimental bounds on the scalar vs. tensor nor- 
malized couplings (upper) and leptoquark exchange helicity 
projection amplitudes (lower panels). T he g rey areas repre- 
sent the information as defined in Ref. [IJ], while the lines 
represent the limits resulting from the present experiment. 
Decreasing line thickness as well as intensity of the grey areas 
correspond to 1-, 2- and 3-sigma confidence levels. 

Similar constraints can be imposed on the selectron ex- 
change couplings (Alii, Aiii) in the minimal supersym- 
metric standard model (MSSM) with R-parity violation. 
Adopting conventions used in Ref. [ssj . the amplitude of 
the selectron (e) exchange between the quark and lepton 
can be written as: 



selectron exchange to the scalar coupling of the beta- 
decay and to the R and N correlations can be written 
as: 



AiiiA-ii 



i? = -A 
A^= -A 



Ami 

1=2,3 <^i'i 



2V2ImCs 



Gj.K,d5y(l + 3A2)' 

2V2ReCs 
GFK,d.9y(l + 3A2)' 



(38) 

(39) 
(40) 



where gs is given by the neutron, proton and light quarks 
masses: 



Mn - Mp 
ruu - ma 



0.49 ±0.17, 



(41) 



Vud = 0.97425(22) is the CKM matrix element, A = 
-1.2694(28) and Gf = 1.166364(5) x 10"^ GcV-^ H. 

The presently best direct constraint for the imaginary 
part of the scalar interaction obtained in the present ex- 
periment improves significantly the limits on the combi- 
nation of coupling constants leading to the updated ver- 
sion of the exclusion plot presented in Ref. [53] (Fig. [331) • 

The real part of this combination was also accessible 
in this experiment via the measurement of the N correla- 
tion. The achieved accuracy, even if slightly better than 
that for _R, can not compete with the much more pre- 
cise data adopted from the compilation of superallowed 
Fermi nuclear beta decays. 



FIG. 33: Experimental bounds on the real vs. imaginary 
combined couplings of MSSM with R-parity violation. The 
grey areas represent the information as defined in Ref. [2l| ]. 
while the lines represent the limits resulting from the present 
experiment. Decreasing line thickness and intensity of the 
grey areas correspond to 1-, 2- and 3- sigma confidence levels. 



yMe, 



E 

1=2,3 



Alii A' 
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u(l-H75)d -6(1-75)^^6, (37) 



where m^^ is the slepton mass assumed to be equal to 

, the contribution of the 



100 GeV. According to Ref. 



VI. 



CONCLUSION 



The measurement of the transverse polarization com- 
ponents of electrons from the decay of free polarized neu- 
trons has been carried out successfully. This was the 
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first experimental determination of the R correlation co- 
efficient in neutron /3-decay and, to our knowledge, also 
the first observation of the finite value of the N correla- 
tion, an effect of the final state interaction in the neutron 
/3-decay. The obtained results allowed a significant im- 
provement in constraining the relative strength of exotic, 
scalar-type weak interaction and related parameters in 
standard model extensions with leptoquark exchange and 
in the MSSM with R-parity violation beyond the limits 
from all previous measurements. 

The most important feature of the experimental setup 
which made this possible was the ability to fully recon- 
struct momenta of low energy electrons before and after 
the backward Mott scattering which served as the elec- 
tron polarization analyzer. 

Further, even substantial improvement in the statisti- 
cal accuracy of the determination of R and A'' correlation 
coefficients can be achieved in an experiment based on 



this principle, provided that a substantial increase of the 
solid angle acceptance is attained. 
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